Logarithmic corrections to finite size spectrum of SU(N) symmetric quantum chains 
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We consider SU(N) symmetric one dimensional quantum chains at finite temperature. For such 
systems the correlation lengths, ground state energy, and excited state energies are investigated 
in the framework of conformal field theory. The possibility of different types of excited states are 
discussed. Logarithmic corrections to the ground state energy and different types of excited states in 
the presence of a marginal opeartor, are calculated. Known results for SU (2) and SU (4) symmetric 
systems follow from our general formula. 
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One of the fundamental models of solid state physics is 
the Heisenberg model of insulating magnets. In the one- 
dimensional case ("spin chains"), the spin- 1/2 Heisen- 
berg models have been studied extensively: most of 
our understanding of their quantum critical behavior is 
based on the Bethe- Ansaiz. solution for the ground state 
and excitation spectrum,Ero mapping to the Sine-Gordon 
theory,S non-abelian bosonization,cl and mapping to the 
sigma model.tj Although spin- 1/2 Heisenberg chains are 
SU(2) symmetric systems, fruitful generalizations have 
been accomplished in two different directions: (a) enlarg- 
ing the representation of SU (2) group to study quantum 
chains with higher spins, and (b) introducing higher sym- 
metry groups such as SU(N). 

Here we consider generalizations of the type (b) and in- 
vestigate how higher symmetry affects the ground state 
properties [Eqs. and finite size spectrum of 

quantum "spin" chains. Earlier studies of Afflecktl show 
that any one dimensional system with SU (JV) symmetry 
is critical, and at the very low energy scale these mod- 
els are equivalent to (N — 1) free massless bosons. This 
system of free bosons - when viewed in the framework 
of two dimensional conformal field theory, are the pri- 
mary fields of the SU(N) k =i WZNW model. Adopting 
this model we give an explicit derivation of logarithmic 
corrections to finite size spectrum of SU(N) symmetric 
quantum chains. Logarithmic shifts in excited states en- 
ergy levels have been theoretically observed for SU(2) 
and SU(4) symmetric systems away from the T = 
quantum critical point .0il3 Known results for N = 2 
and N = 4 follow from our general formulas [obtained in 
Eqs. @-(§f) and the paragraphs below Eq. Q]. 

An one-dimensional SU(N) symmetric quantum chain 
of length L is described by the HamiltonianEj 
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(do being the lattice spacing) and S A are the (N 2 — 1) 
generators of the SU (N) Lie algebra at each lattice site j. 
For convenience, the interaction strength and the lattice 
spacing clq have been set equal to one in Eq. (|l|). At 
each site j, the generators Sf can be represented by N 
"flavors" of fermions, ip a j [a = 1, N] 
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S f=J2 ^(T A ) b a i, bj -I/N, 
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where I is the identity operator, and T s are a com- 
plete set of (N 2 — 1) traceless normalized matrices so 
that Tr[T A T B ] = (1/2)S AB . Eq. (§) satisfies the con- 
straint that at each site the total number of fermions is 
conserved, i.e. Y^a=i^j = 1 - 

For such a fermionic system the theory can be 
bosonized using non-abelian bosonization at low tem- 
peratures. n In the continuum limit, the bosonized 
Hamiltoniano (H c g) can then be written in terms of the 
Kac-Moody currents, 
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Jd x [j£j£ + J A J A + 2J A J A ], 
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where the normal ordered Kac-Moody currents for the 
left (with Fermi momentum kp < 0) and the right mov- 
ing (with kp > 0) fermions are defined as, 
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Here n = L/clq is the total number of discrete points 



At zero temperature (T = 0) the interaction term, 

^^J^ 1 J A Jr , m Eq. (||) renormalizes to zero, and the 
sum of first two terms in the Hamiltonian that are 
quadratic in left and right moving currents, corresponds 
to the SU{N) k=1 WZNW model. The fundamental uni- 
tary N x N matrix field g of the WZNW model is given 
by, 
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9b = (const) : i/j^ipm 
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The field g transforms as the fundamental representation 
of SU(N) L x SU (N) R which describes the exact symme- 
try of the Hamiltonian in Eq. (^) at zero temperature. It 
is known that this fermionic theory is equivalent to a the- 
ory of (N — 1) free massless bosons at the criticality with 
velocities v s ; they correspond to (N— 1) excitation modes 
of the SU(N) symmetric rauantum chain that oscillates 
at different values of kpB Furthermore, these oscillat- 
ing modes are primary fields of the SU(N)k=i WZNW 
modeLand their scaling dimensions (A p ) can be obtained 
frorrO 
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For SU(N), C p (p = l, 2, • • • , N - 1) is the eigenvalue of 
the Casimir operator in the p-th fundamental represen- 
tation [having a Young tableau with p boxes in a single 
column], and C ac ij = N is the eigenvalue of the Casimir 
operator in the adjoint representation [having a 2-column 
Young tableau with (N — 1) boxes in the first column 
and 1 box in the second column]. In a highest weight 
(A) representation of SU(N), the corresponding Casimir 
eigenvalue is given bycJ 



C A = 
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where 9 is the highest weight corresponding to the ad- 
joint representation and is normalized to 1, m is the 
total number of boxes in the Young tableau with r 
rows of length b\, 62, ■ • ■ , W a and c columns of length 
ffli, 02, • ■ • , a Co . Using this formula we find C p = p(N — 
p)(N + l)/2N and C a( jj = N. Hence, the scaling dimen- 
sions (A p ) of primary fields of the SU(N) k =x WZNW 
model are given by 
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For example, in the case of SU(4) the three oscillating 
components have scaling dimensions (3/4, 1, 3/4) [for 
p = 1, 2, 3]. For SU{N), the mode (dominant) that 
oscillates at hp = 2tt/N has a scaling dimension (1—1/N) 
[p = 1 or (N — 1) in this case]. 

Finite-size corrections to Heisenberg chain with SU (2) 
and SU(4) symmetry has been studied using confor- 
mal field theory.QO The relevance of studying finite size 
chains are twofold. One can not only compare the the- 
oretical results with numerical simulations and experi- 
ments which are limited to finite size of the system but 
also can study the finite temperature behavior of the sys- 
tem by identifying the finite size in the imaginary time 



direction, that corresponds to finite temperature. To ob- 
tain the finite-size corrections of a ID chain of length 
L and with periodic boundary conditions we first intro- 
duce a conformal mapping from the infinite plane (with 
coordinate z) to the cylinder (with coordinate w) via 
w = (L/2ir)hiz. Identifying the length to be inverse 
of thep-temperature [L = v s /T] the finite temperature 
resultaB of the ground state energy Eq can be general- 
ized to the SU (N) symmetric system, 



E o (T)=E o (0) 



tvT(N-I) 
6v s 



(9) 



Here Eq(0) refers to the ground state energy at zero 
temperature. The thermodynamic quantities like spe- 
cific heat and entropy can now be obtained by taking the 
appropriate derivatives with respect to the temperature. 

Other quantities of interest are the finite temperature 
corrections to the correlation lengths (£) of the different 
modes. These inverse of the correlation lengths, £ -1 are 
signature of energy gaps, (E n — Eq) between the ground 
state and the lowest lying excited states (E n ) that are 
created by finite temperature of the system. Using the 
general formula for the scaling dimension [Eq. (||)], we 
obtain t;^ 1 of the p-th staggered mode, 
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The temperature dependence of the correlation lengths 
is in fact modified by logarithmic corrections in the pres- 
ence of marginal operators in the theoryO The generic 
form of the Hamiltonian at the critical point containing 
a marginal operator cf>(x,t) is 



H = H* +g Q / dx 



(11) 



where g Q is the coupling constant and H* is the Hamil- 
tonian at the fixed point. In our case [Eq. (Q)], the nor- 
malized marginally irrelevant operator is 
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For such a marginally irrelevant operator, the Hamilto- 
nian at the critical point (T = 0) becomes equal to the 
fixed point Hamiltonian. In Eq. (I4), D is the normaliza- 
tion constant to be determined from the two-point cor- 
relator of (f>. The operator product expansion (OPE) of 
the J^tR) w ith any normalized Kac-Moody primary field 
X is given by0 
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where the operators Jg ^ and Jqr are the genera- 
tors of the global SU(N) L x SU(N) R transformations, 
and satisfy the characteristic equations, Jq L \x( z ')) — 
-T£\ X (z')) and j£ R \ X {z')) = \x{z'))Tl ' Note that 
J2a(Jol) 2 an d 12a(Jor) 2 are ^ ne Casimir operators of 
SU(N) L and SU(N)r groups, respectively. The two- 
point correlators of the left currents (for k = 1) are, 



{^(z)JE(z')) 



Tr[TfTf] 
"47T 2 (z- z') 2 
S AB 
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Using these results we explicitly calculate: 
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and then compare it to the standard conformal field the- 
ory result i.e. ((j>(z, z)4>{z' , z')) = \z — z'\~ 2 \z — z'\~ 



obtain the value of the constant, 
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Thus the normalized irrelevant marginal operator is given 
by, 
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Perturbation to the normalized excited state (</)„) ene& 
gies due to the marginal operator can now be calculatedtll 
from 



5(E n - E ) = g a Jdx (<p n \4>\4> n ), 



(20) 



where <p and <p n are Virasoro primary fields generated by 
applying Fourier modes of J A and J R on Kac-Moody pri- 
mary fields. For large length (equivalently, small temper- 
ature), we may replace the coupling go by its renormaL 
ization group improved value (upto the log — log term) ,t3 
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21n(T /T) 



ln[ln(T /T)] 
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Here To is the model dependent parameter of the system 
and the coefficient b is defined via the following 3-point 
correlator, 



Zl)4>(z2,Z 2 )<p(z 3 ,Z 3 )) = -~b/\z 12 \ 2 \z23\ 2 \z 13 \ 2 . 

Substituting this in Eq. (|2(i| ) we obtain, 
2?rT 



(22) 
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v s ln(T /T) 
1 

1 



21n(T /T) 



ln[ln(T /T)] 
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The coefficient b n is again defined through the 3-point 
correlator, 



<f> n (zi, Zi)4>(z 2 , z 2 )cj> n (z 3 , z 3 )} 
-b n /\z 12 \ 2 \z 23 \ 2 \z 13 \ 2 ^- 2 . 



(24) 



Here x n is the scaling dimension of the Virasoro primary 
field (j> n . Substituting Eq. ( |l9| ) in Eq. ( p4| ) and using 
the OPE's as in Eqs. (|l3|) and (0) it follows that b n is 
directly proportional to the sum of the product of the 
eigenvalues of the generators J A L and Jq R , 
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To evaluate J2a t l t r> 

we observe that the full symme- 
try, SU(N) l x SU(N)r, of the quantum chain at the crit- 
ical point is broken by the presence of the marginal op- 
erator, (j)(z,z). Only the diagonal SU(N) C SU(N) L x 
SU(N) R is an exact symmetry of the quantum chain. 
Under this subgroup, the representation Vl <£> Vr of 
SU{N)l x SU(N)n decomposes into direct sum of var- 
ious irreducible subrepresentations. If an excited state 
(l^n)) belongs to a highest weight subrepresentation, 
V C Vl ® Vr and C is the corresponding Casimir in- 
variant of the diagonal SU (N) in V, then we haveO 



TV -1 

E t£t a = -[c-c l -Cr] 

A=l 1 



(26) 



where Cl and Cr are the Casimir invariants of SU(N)l 
and SU(N)r in the highest weight representations Vl 
and Vr, respectively. Therefore, using Eqs. ( p5| ) and ( |2^ ) 
we find 
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[C-C L - Cr] . 



(27) 



The above formula may also be used to evaluate the 
renormalization group coefficient, b [Eq. (B2)]: since 
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<f>(z, z) is a Virasoro primary field of conformal dimen- 
sions (1,1), we set (f) n (z,z) = </>(z,z) and x n = 2 in 
Eq. (Hj), and hence, b n = b. This can be seen as fol- 
lows. The Virasoro primary fields j£ and of confor- 
mal dimensions (1,0) and (0, 1) transform as the adjoint 
representations, and of SU(N), and since V R * 
is conjugate to V£ \ the direct sum decomposition of 

V L di ® V B di under the diagonal SU(N) C SU(N) L x 
SU(N)n must contain a unique singlet. Hence, the Vi- 
rasoro primary field <p(z, z) in Eq. (19) transforms as this 
singlet representation and we have C — 0, Cl = Cr = N 
in Eq. ([27|). This implies 
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For example in the case of Heisenberg spin-chain with 
SU(2) symmetry, b = 4/\/3 and for the spin-orbital 
model with SU(4), b = 8/VTE. Our result for SU(4) is 
new. Together with Eq. (pij) , the constant b also deter- 
mines the correction of 0(g%) in the ground state energy, 

E Q (T) - E (0) = [(N - 1) + 2n 3 bg 3 ] . (29) 



To determine the logarithmic shifts in the excited 
states energy levels, we need the ratio 2b n /b in Eq. (23). 
From Eqs. (p7|) and (E8h we get, 
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= ~[C-C L -C R ] 



(30) 



To evaluate b n (and hence 2b n /b) we must know the 
excited states. In SU(N) invariant quantum chains, the 
low lying excited states (\4> n }) correspond to (V— 1) pri- 
mary fields of the SU(N) k=1 WZNW model with the 
scaling dimensions, A p . These fields transform as (q, q) 
representations of SU(N)l x SU(N)r, where 



Q = 



N(N -!)■■■ (N -p+1) 



p\ 



is the dimension of the p-th fundamental representation 
of SU(N) for p = 1, 2, • • • , N - 1. For instance, the 
lowest excited states correspond to the fundamental pri- 
mary field g with the scaling dimensions, Ai = (1 — 1/JV). 
This field transforms under the (V, N) representation 
which decomposes under the diagonal SU (N) into the ad- 
joint and singlet representations. For the excited states, 
Tr [gT' 4 ], belonging to the adjoint representation, we 
have C = N, C L = C R = (N 2 - 1)/2N which implies 
b n = -1/(NVN 2 - 1) and 2b n /b = -1/N 2 . For exam- 
ple, in the case of SU(2), this ratio 2b n /b = — 1/4 J3 and 
for SU(3) and SU(4) they are -1/9, and -1/16 respec- 
tively. 



For the excited state, Tr g, belonging to the singlet 
representation we have C = 0, C L = C R = (N 2 - 1)/2N. 
In this case, b n = V N 2 -1/N and the 2b n /b = 1 - 
1/N 2 . This result is new. In case of SU(2), J t,he value 
2b n /b — 3/4 has been previously obtainedQ but for 
SU(3), 2b n /b = 8/9 and for SU(4), 2b n /b = 15/16 are 
predictions from our general formula. 

We consider one more application of the formula ( |30| ) 
of current interest - the SU(4) symmetric quantum chain 
described by the SU(4)k = i WZNW model. In this case, 
to compute logarithmic corrections to the excited states 
energy we note that there are three primary fields with 
scaling dimensions A p = 3/4, 1, 3/4 for p = 1, 2, 3 
respectively, as seen from Eq. (0). The case of p = 1 
[and p — 3], as discussed above, is the fundamental field 
g [and its hermitian conjugate g] which transforms under 
the (4, 4) [and (4, 4)] representation of S77(4) L x SU(4) R . 
From Eq. jl0|), the next lowest energy excited states cor- 
respond to the primary field operator (denoted by \I/) 
with A2 = 1. The field <J> transforms under the (6, 6) 
representation of SU (4) l x SU (4) R , This (6, 6) represen- 
tation decomposes as direct sum of a singlet, an adjoint 
and a 20-dimensional representation (as in Fig. [l]) under 
the diagonal SU{4). 
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FIG. 1. Young tableau for the decomposition of (6, 6) rep- 
resentation of SU(4:). The number in the parenthesis denotes 
the dimension of the corresponding representation. 

We now compute the ratio 2b n /b for the excited states 
corresponding to 20-dimensional representation which 
has a Young tableau with 2 rows and 2 columns. For 
this representation, the Casimir invariant C in Eq. ( |30| ) 
is obtained from the formula (Q): we find C = 6, and 
Cl = Cr = C p= 2 = 5/2. Thus, for the excited states 
corresponding to in the 20-dimensional subrepresenta- 
tion of (6, 6), we have 2b n /b = —1/4. 

In summary, we have studied finite size spectrum for 
one dimensional SU(N) symmetric quantum chains us- 
ing both conformal field theory and representation theory 
of SU (N) . We have calculated in general the scaling di- 
mensions of all the oscillating modes, and obtained the 
ground state energy as well as correlation lengths of the 
staggaered modes for a finite size system with SU(N) 
symmetry. Possibilities of different types of excited states 
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are also briefly discussed and a general formula to com- 
pute the logarithmic correction to the excited state ener- 
gies has been derived. The existing results for N = 2, 4 
agree with the predictions from our general formula. 
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